THE ONSET OF CHAOS IN PULSATING VARIABLE STARS 
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ABSTRACT 

Random changes in pulsation period occur in cool pulsating Mira variables, Type A, B, and C semiregular 
variables, RV Tauri variables, and in most classical Cepheids. The physical processes responsible for such 
fluctuations are uncertain, but presumably originate in temporal modifications of the envelope convection in 
such stars. Such fluctuations are seemingly random over a few pulsation cycles of the stars, but are dominated 
by the regularity of the primary pulsation over the long term. The magnitude of stochasticity in pulsating 
stars appears to be linked directly to their dimensions, although not in simple fashion. It is relatively larger 
in M supergiants, for example, than in short-period Cepheids, but is common enough that it can be detected 
in visual observations for many types of pulsating stars. Although chaos was discovered in such stars 80 years 
ago, detection of its general presence in the group has only been possible in recent studies. 
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A well known problem in variable star studies is that 
it is impossible to predict exact moments for light maxi- 
mum in some late- type pulsating variables, such as Miras 
and semi-regular variables, or to predict their amplitude 
on any given cycle (see Fig. [T] for AAVSO (American 
Association of Variable Star Observers) observations of 
the name star, o Ceti). The cyclical light patterns dis- 
played in such stars are reasonably well defined over long 
time intervals and can be approximated closely with lin- 
ear ephemerides, but the regularity of their pulsation is 
typically marked by other effects best revealed through 
careful O-C analysis. 

A common complication is that of "random" fluctu- 
ations in pulsation period for a star from one cycle to 
another. Many years ago Eddington & Plakidis (1929) 
developed an interesting technique for establishing the 
importance of random fluctuations in pulsation period 
for Mira variables, and it has been revived frequently in 
recent years (Percy & Hale 1998; Percy & Colivas 1999, 
Percy et al. 1993, 2003, 2007) in order to establish the 
importance of random changes in period for other Mira 
variables as well as for other types of cool and hot pul- 
sating variables. 

The technique has been described previously in these 
pages (Turner & Berdnikov 2001), and involves comput- 
ing, without regard to sign, the average accumulated 
time delays (u(x)) between light maxima separated by 



x cycles. If the deviations in the observed times of light 
maxima from their predicted times are dominated by 
random fluctuations in period, then the data for all avail- 
able observed light maxima should display a trend de- 
scribed by: 

{u(x)) 2 = 2a 2 +xe 2 , 

where a is the average uncertainty in days for estab- 
lished times of light maxima and e is the magnitude of 
any random fluctuations in period. The technique could 
alternatively be formulated to represent e in terms of 
phase offset. A schematic representation of expecta- 
tions for a "typical" pulsating variable is presented in 
Fig- 12 But such expectations are never matched exactly 
in practice, since the dominant pulsation in such stars 
forces the random factor in the pulsations back into a 
regular pattern of variability after ~ 50 — 200 cycles. 



2. Different Types of Pulsating Stars 

Although the Eddington-Plakidis test was devised 
specifically to analyze random fluctuations in pulsa- 
tion period for Mira variables, and later applications by 
Percy et al. (1993), Percy & Hale (1998), and Percy & 
Colivas (1999) included "Mira-like" stars (Miras, Type 
A and B semiregulars, and RV Turi variables), the same 
technique should be valid for all types of stars. The 
technique was also applied to short period pulsators of 
both Population types by Percy et al. (2003, 2007), was 
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Fig. 1. — The visual light curve of Mira between JD 2435000 and JD 2440000 from AAVSO measurements. 
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Fig. 2. — A schematic representation of an Eddington- 
Plakidis test for a pulsating variable with a randomness fac- 
tor of e = 2.7 days and uncertainties of a = 0.15 day in 
measured times of light maximum. 



extended to Cepheids in a number of studies (Turner & 
Berdnikov 2001, 2004; Abdel-Sabour Abdel-Latif 2004; 
Berdnikov et al. 2004, 2007, 2009a, 2009b; Berdnikov 
2010, Berdnikov & Stevens 2010). Mostly negative re- 
sults were found for short period Cepheids, but that is 
because individual light maxima are rarely observed for 
such stars and the individual times of light maximum 
tabulated in O-C analyses usually refer to data obtained 
over many adjacent cycles about the one cited. Since 
random fluctuations in period exist over several pulsa- 
tion cycles, their effects on times of light maximum can 
easily be confused with other sources of scatter in the 
light curves when the data are averaged over many cy- 
cles. 

Confirmation of that conclusion has recently come 
from space observations of Cepheids (Berdnikov 2010; 
Berdnikov & Stevens 2010). Since weather problems and 
limited observing windows are generally not a problem 
for space observations, it is possible to observe many 
consecutive light maxima for short period Cepheids and 
to derive observed times of light maxima from appli- 



cations of Hertzsprung's method. The results indicate 
that random fluctuations in period are relatively com- 
mon even for short period Cepheids, and most likely 
apply to all radially pulsating stars, at least to some 
extent. 

The same feature also appears in at least one SRC 
variable, the pulsating M3 la supergiant BC Cyg, a 
star for which the observed times of light maximum can 
vary by ±84 days from those predicted by a quadratic 
ephemeris accounting for its long-term period decrease 
(Turner et al. 2009). An Eddington-Plakidis test for the 
star is replicated in Fig.[3j which demonstrates that ran- 
dom fluctuations in pulsation period for the star domi- 
nate over small cycle counts, whereas the dominant reg- 
ular pulsation appears to return the variability to a regu- 
lar pattern after about 12 cycles or so (~ 23 years in the 
present case). In pulsating stars of shorter cycle length 
such regularity returns much sooner, a matter of a few 
years in the case of short period Cepheids (Berdnikov 
2010). 

Evidence for random fluctuations in pulsation period 
on short time scales can also be seen in the results of 
Poleski (2008) for several Cepheids in the Large Magel- 
lanic Cloud. The evidence is revealed by offsets in the 
observed times of light maximum from cycle to cycle 
in Poleski's O-C diagrams, although the deviations are 
typically rather small in comparison with much larger 
deviations observed over longer time intervals, where the 
evolutionary changes in mean radius become dominant 
(Turner et al. 2006). 

3. Parameterizing the Randomness Factor 

In their original tests on o Ceti (Mira) and x Cygni, 
Eddington & Plakidis (1929) noted that the observed 
random fluctuations in period for both stars amounted 
to about the same amount, 1.35% to 1.39% of the pul- 
sation period. They clearly understood the importance 
of the star's pulsation period to the magnitude of the 
stochastic processes producing the random fluctuations 
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Fig. 3. — An Eddington-Plakidis test for the M3 la super- 
giant variable BC Cyg. 
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Fig. 4. — The observed trend of increasing randomness fac- 
tor "e" with increasing pulsation period P. 



in period. Yet the cycle length for any pulsating star also 
depends directly on stellar radius through the period- 
radius relations applying to every type of pulsating star. 
Therefore, a better parameter for describing the stochas- 
tic processes arising in pulsating stars is the ratio e/P, 
although that is only to first order. The parameter e/P 
must be independent of radius if the variables obey sim- 
ilar period-radius relations. 

Results to date for all Eddington-Plakidis analyses of 
pulsating stars (updated from Turner et al. 2009) are 
illustrated in Figs. 0] and [5l and confirm such an assump- 
tion. Note that Fig. 2] contains additional data from 
earlier versions plotted by Turner & Berdnikov (2001) 
and Turner et al. (2009). As seen in Fig. [5j the param- 
eter e/P is indeed relatively independent of pulsation 
period, in other words independent of stellar radius, al- 
though there may be an additional trend with period, 
perhaps reflecting the increasing dominance of envelope 
convection with decreasing stellar surface temperature 
in such stars. The parameter e/P has a mean value of 
0.0136±0.0005 (±0.0069 s.d.), which matches the results 
of Eddington & Plakidis (1929) more than 80 years ago. 

The actual variation of randomness factor e with pul- 
sation period P may be somewhat more complicated. 
For example, most of the stars with the smallest val- 
ues of e/P are short period Cepheids of relatively high 
mean effective temperature, whereas the star of longest 
pulsation period, BC Cyg, may have a larger value of 
e/P than indicated in Figs. [4] and [5j if one considers 
the variation of (u(x)) 2 over small cycle differences x in 
Fig. El The M3 la supergiant BC Cyg is, of course, the 
star of lowest mean effective temperature in the sample. 

As noted earlier, the actual observed trends in the 
computed values of (u(x)) 2 for all pulsating stars tested 
to date initially increase directly in proportion to in- 
creasing differences in cycle count predicted (Ed- 
dington & Plakidis 1929, see also Fig. [5]). At larger cy- 
cle differences, however, the trend reverses as the dom- 
inant pulsation in such stars reimposes its regularity in 



the observed times of light maximum. Observed values 
of (u(x)) 2 gradually become much smaller for large cy- 
cle differences, a characteristic also noted by Eddington 
and Plakidis in their original study of two Mira vari- 
ables and seen in almost all of the Eddington-Plakidis 
tests cited earlier. The stochastic fluctuations in period 
that appear as a common feature in the cycle lengths 
of nearly all pulsating stars are therefore dominated by 
the regular pulsation in such stars. As noted above, 
the physical processes responsible for such characteris- 
tics are uncertain, but presumably originate in temporal 
modifications of envelope convection in such stars. 

4. Discussion 

The fact that random fluctuations in period are ubiq- 
uitous for all pulsating variables has important conse- 
quences. Standard O-C analyses of Cepheids, for ex- 
ample, will always display scatter in the individual O-C 
data, unless they are averaged over many adjacent cy- 
cles. Even then, the non-photometric source of scatter 
in light curve data points for individual cycles must also 
introduce a small source of uncertainty in the resulting 
O-C data based on light curves constructed from ob- 
servations averaged over the same cycles. Likewise, the 
detection of sizable random fluctuations in period for 
long period variables means that the Predicted Dates of 
Maxima and Minima of Long Period Variables issued 
regularly by the AAVSO must necessarily be inexact. 
Fortunately such predictions are generally issued within 
a cycle or two of the predicted dates, so they are likely to 
be only a few days off because of the stochastic processes 
occurring in the envelopes of such stars. 

The AAVSO collection of data for Miras and long pe- 
riod variables was the basis for studies of random fluctu- 
ations in their pulsation periods by Percy et al. (1993), 
Percy & Hale (1998), and Percy & Colivas (1999), so it 
seems clear that such changes can be detected in long 
period variables from simple eye estimates. An interest- 
ing question to ask is whether or not such effects can also 
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Fig. 5. — The nearly negligible trend of the parameter e/P 
with pulsation period. 



be detected from simple eye estimates for Cepheid vari- 
ables. The lead author recently had an opportunity to 
address that question while teaching an undergraduate 
course in astronomy. A decade previously, Turner (1999, 
2000) demonstrated a simple procedure for obtaining 
precise estimates of magnitude for bright Cepheids from 
unaided observations by eye. The original observations 
from 1998-99 have been used frequently since then for 
instructional purposes, and the procedure was revived 
in October 2009 to provide a reference data set for a 
student attempting to follow the technique on his own. 
The results are illustrated in Fig. [6j where the data are 
phased using the same ephemeris adopted in 1998-99. 

The observed light maxima for S Cephei in winter 
1998-99 and fall 1999 fell very close to the times pre- 
dicted from the adopted ephemeris for the star, and 
the same conclusion applies to the observations from 
2009. The Cepheid has a well-established period de- 
crease (Turner 1999, 2000; Turner et al. 2006), so the 
times of light maximum should occur slightly earlier in 
2009 than they did in 1998-99. But the O-C offset is ef- 
fectively nil between the two dates, and no evidence for 
chaotic effects in pulsation period can be distinguished 
from our eye estimates for the star. The results of Berd- 
nikov (2010) and Berdnikov & Stevens (2010) for other 
short-period Cepheids are consistent with such conclu- 
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Fig. 6. — Observations by Turner with unaided eye of 8 
Cephei from 1998 to 2009. 
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